Evidence for effective thermal boundary resistance from magnon/phonon disequilibrium 
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We use the time-resolved magneto-optical Kerr effect (TRMOKE) to measure the local temperature and heat 
flow dynamics in ferromagnetic SrRuOs thin films. After heating by a pump pulse, the film temperature decays 
exponentially, indicating that the heat flow out of the film is limited by the film/substrate interface. We show 
that this behavior is consistent with an effective boundary resistance resulting from disequilibrium between the 
spin and phonon temperatures in the film. 
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Propagation of heat in magnetic nanostructures is a sub- 
ject of importance for both fundamental and practical reasons. 
Regarding applications, devices based on spin-current torque 
dissipate significant amounts of energy in the process of flip- 
ping spins [1,2], and successful operation requires that energy 
escape from the structure in the form of heat. From a funda- 
mental point of view, the factors that limit the flow of heat 
out of magnetic nanostructures are not well understood. For 
example, macroscopic (Fourier) modeling of heat flow is not 
applicable when sample dimensions are smaller than phonon 
or other quasiparticle mean-free-paths. Heat propagation is 
more complex in ferromagnetic nanostructures because ther- 
mal energy is shared over spin, electron, and phonon degrees 
of freedom (DOF). The net rate of heat flow depends on the 
thermal diffusivity of each DOF, the rate at which energy is 
shared between them, and finally, the extent to which each 
DOF is confined to the nano structure. 

A deeper understanding of heat flow in magnetic nanos- 
tructures requires experiments that explore the dynamics of 
local temperature change over a wide range of parameters, for 
example, bath temperature and sample dimensions. Optical 
techniques such as as thermal reflectometry have proven use- 
ful in measuring the local temperature in nanostructures. Dy- 
namics are probed effectively in time-resolved versions of this 
technique, in which a focused laser pulse heats the sample and 
a second time-delayed pulse probes the temperature via the 
associated variation in optical reflectivity One draw- 

back of time-resolved thermal reflectometry (TRR) is that it 
is limited to structures whose dimension, L, is greater than 
the penetration depth, Ap, of the laser light. If L <C Ap most 
of the pump light is absorbed in surrounding materials rather 
than nano structure, which leads to changes in reflectivity not 
associated with heating of the sample. 

When the structures of interest are magnetic, the time- 
resolved magneto-optical Kerr effect (TRMOKE) offers po- 
tential advantages over TRR. This measurement uses the tem- 
perature dependence of the Kerr rotation, A6k, rather than 
reflectivity, as the thermometer. The most significant advan- 
tage relative to TRR is that A6k arises only from the sam- 
ple under study, even when L <C Ap. A second advantage is 
that TRMOKE can be a more sensitive probe of temperature 



in magnets in which strong spin-orbit coupling gives rise to 
large Kerr effects. 

In this work we use TRMOKE to study the dynamics of 
heat transport in ferromagnetic thin films of SrRuOs (SRO). 
SRO is ideal for studying heat transport in nanoscale ferro- 
magnets for several reasons. First, high-quality, epitaxial films 
with thicknesses ranging from 2-200 nm can be grown on 
SrTiOs substrates by pulsed-laser deposition. Moreover, SRO 
has an unusually large Kerr coefficient that allows measure- 
ment of local changes in temperature of a few Kelvins in films 
that are only 2 nm thick [7] . Finally, as we discuss further be- 
low, the relevance of the magnetic DOF in thermal transport 
is accentuated in SRO compared with the more widely studied 
3d ferromagnets such as Fe, Ni, and Co. 

Fig. la shows the TRMOKE signal observed in an SRO 
film of thickness 200 nm, at several temperatures, T, be- 
low the Curie temperature of 150 K. The curves show A6x, 
as function of time t, after absorption of the pump pulse. 
ABx(t) is proportional to the z-component of M because 
the probe beam is at near normal incidence. The 200 GHz 
oscillations observed when T < 80 K and t < 10 ps cor- 
respond to damped precession of the magnetization, M 
These oscillations are stimulated by a sudden change in the 
direction of the magnetic easy-axis, caused by absorption of 
the pump beam. The fact that AQxit) approaches a nonzero 
value when the oscillations have died down indicates that M 
has spiralled into alignment with a new local easy-axis direc- 
tion. Below we show that the value that AQk reaches when 
the oscillations have died away provides a calibrated measure 
of the increase in local T induced by absorption of the pump 
pulse. 

In Fig. lb we plot AQk at a fixed time delay of 25 ps as a 
function of T, for three different values of the pump laser in- 
tensity. (Shown in the inset to Fig. lb is the T dependence of 
the equilibrium Kerr rotation, 6^, for the same 200 nm SRO 
sample). The peak in AQ k{T) near the Curie temperature 
has been reported previously in a survey of TRMOKE in a va- 
riety of perovskite and spinel ferromagnets f^. In this work 
we focus on a new feature, which is the increase of A0x ob- 
served as T ^ 0. 

If AQk for t > 10 ps is caused by local heating, we can 
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FIG. 1: (a) Thermally-induced change in Gk for different values of 
pump pulse intensity. Inset: Sk{T) (b) Measured change in tem- 
perature in response to heating pulse. Black lines represent a fit to 
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use &k{T) to calibrate the conversion from the TRMOKE 
signal to the local temperature, Ti{t). To test whether a local 
heating mechanism is consistent with the TRMOKE data, we 
consider the expected scaling of AQ K{t) with bath tempera- 
ture, T, and energy deposited by the laser pulse, ^. After $ 
has been converted to heat, Ti is determined by the relation, 
^ = U{Ti) - U{T), where U{T) is the internal energy. To ex- 
hibit the scaling relation, we express Ti in terms of the inverse 
of U{T), 



Ti{^,T) = U-\^^U{T)). 



(1) 



According to ©, plots of Ti{^, T) at fixed T will fall on a 
single curve, /7~^(<l>), if shifted along the abscissa by U (T). 

In Fig. 2a (left panel), we plot T/($, T) at a delay of 25 
ps as a function of <l> for different values T. The values of 
Ti were determined by comparing AQ K{t = 25ps), with the 
calibration data, Qk{T). The result of shifting the curves to 
achieve the best collapse to a single curve is shown in the right 
panel of Fig. 2a. The data are seen to exhibit the scaling prop- 
erty expected for the local heating model. A powerful feature 
of this scaling approach is that we can determine the inter- 
nal energy function U (T) from the magnitudes of the shifts 
required for the data collapse (shown as horizontal arrows in 
Fig. 2a). The internal energy as a function of T determined by 
this procedure is shown in Fig. 2b. When plotted as a function 
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FIG. 2: (a) Unshifted measurements of final temperature as a func- 
tion of pump laser intensity. Lines represent different cryostat tem- 
peratures, (b) Data from (a) shifted along the [/-axis to lie along a 
common curve, (c) Internal energy vs. T. Data show a T^^^ depen- 
dence below 75 K 



of T^/^ the internal energy appears as a straight line with zero 
intercept for T < 75 K. 

The observation that U (T) is proportional to T^/^ (cor- 
responding to specific heat C ex T^/^) indicates that in 
SRO thermal energy is stored predominantly in magnons for 
T < 75 K. This contrasts with most ferromagnets, in which 
Cp » Cm, where Cp and Cm are the phonon and magnon 
contributions, respectively. However, it is a natural conse- 
quence of the reversal in the ordering of the Curie and De- 
bye temperature (Tq) in the two classes of materials, that is 
Tc < Tq in SRO whereas Tc > Tq in the elemental 3d- 
transition metal ferromagnets. 

From the preceding analysis, we conclude that the measure- 
ment of AOk (T, t) provides a direct probe of Ti (T, t) in SRO 
thin films. In the following we explore the dynamics of Ti as a 
function of time delay, substrate temperature, and film thick- 
ness. In Fig. 3a we plot AQxit) for 10 ps < t <500 ps, 
measured at 10 K, for films of different thickness in the range 
from 10-200 nm. The laser power is lowered so that AQk 
is directly proportional to ATi{t), the difference between the 
instantaneous local temperature at the substrate temperature. 
We find that ATi{t) decays exponentially, with a time con- 
stant, r, that increases with sample thickness L. Fig. 3b is 
a plot of r as a function of L, showing that for L > 25 nm. 
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FIG. 3: (a) Semi-log plot of AG k (t) for 10, 30, 50, 125, and 200 nm 
samples (b) Decay times vs. film thickness at 5 K (c) Decay times 
vs. temperature for 5, 50, and 200 nm samples 



r aL + 70 ps, where a 2.7 ps/nm. Exponential decay 
with a r that increases linearly with thickness indicates that 
the interface between the substrate and the SRO acts as bot- 
tleneck that limits heat flow out of the film. In this case the 
decay of ATi is exponential with r = RCL, where R is the 
thermal boundary resistance (TBR). 

In Fig. 3c we plot r/L, (equal to the EC product) as a func- 
tion of temperature for films of thickness 5, 50, and 200 nm. 
For the thickest film, the RC product is almost independent of 
T from 10-75 K, indicating that R oc 1/C over a broad range 
range of temperature. For the thinner samples RC increases 
with decreasing T below about 50 K. 

The TBR phenomenon is well-known for dielectric and 
metallic films on insulating substrates. In these systems TBR 
is understood to result from the scattering or specular reflec- 
tion of phonons incident on the interface [5, 10]. This mech- 
anism implies that R ex {\t\'^Cp)~^, t is the transmission co- 
efficient for phonons at the interface. For most materials the 
total specific heat C is dominated by phonons, in which case 
Cp cancels from the RC product. Thus r essentially indepen- 
dent of T is predicted, in agreement with many experimental 
observations. 

As we have shown above, SRO exemplifies a system in 
which the specific heat is dominated by magnons rather than 
phonons. In this limit, r(T) should vary as Cm{T)/Cp{T), 
which is approximately proportional to T~^/^. Instead, as 
seen in Fig. 3c, r is essentially independent of T, except for 
the thinnest film. 

The lack of dependence of r on T suggests that the bound- 
ary resistance is not simply proportional to C~^, but is some- 
how related to Cm {T). At first glance, this seems unlikely, be- 
cause magnons are confined to the magnetic film and cannot 
penetrate into the STO substrate. This argument fails to con- 
sider the crucial fact that for heat to leave the SRO film, the en- 



ergy stored in the magnon DOF must flow to the phonon DOF. 
For steady-state heat flows, it has been suggested lllllil] that 
disequilibrium between phonons and a confined DOF creates 
a mechanism for TBR different from the conventional phonon 
impedance mismatch picture. SRO/STO is an ideal testing 
ground this new mechanism because the fraction of thermal 
energy stored in the confined DOF is far larger than in the 
case of the metal film/dielectric interface. 

Analysis of the flow of energy among electrons, phonons, 
and magnons requires a three-temperature model (3TM) [13]. 
The 3TM assumes that each DOF is internally equilibrated so 
a temperature can be associated with each "fluid." We assume 
one-dimensional heat diffusion because the diameter of the 
excited region on the sample is ~ 50 /im, much larger than 
both the optical penetration depth (37 nm) and the film thick- 
ness. Within the 3TM, the differential equation governing the 
electron temperature Tg, for example, is: 



dt 



d 



dz dz 



' Cqp^Tqp - 



- Gp^rn.ATpn 



-S. (2) 



Here z is the normal coordinate, Ke is the electron thermal 
conductivity, and S{z^t) represents heat input from the laser. 
Gep (Gem) govcms the rate of heat flow from the electron to 
the phonon (magnon) DOF and AT^p {AT em) is the electron- 
phonon (magnon) temperature difference. Permutation of the 
subscripts generates equations for the time-evolution of the 
phonon and spin temperatures, Tp and Tm, respectively. 

To develop a picture of the heat flow dynamics, we have 
used a Dufort-Frankel finite-difference scheme [6] to integrate 
the 3TM equations. In the course of analyzing the 3TM equa- 
tions, we have found that a 2TM model is sufficient to capture 
the essential physics. The electron and magnon DOF can be 
lumped together as single confined DOF, with thermal con- 
ductivity Kc = Km ^ i^e^ spccific heat Cc = Cm + Cg, and 
confined fluid temperature Tg. We note for SRO in the range 
10 K< T <70 K that i^c - and Cc ^ Cm- For afl calcu- 
lations we have assumed that phonon thermal conductivity in 
the substrate is much larger than in the film, which enforces 
ATi at the film/substrate interface. Furthermore, we as- 
sume no heat flow from the surface of the film to vacuum, 
corresponding the boundary condition that AT"{z) vanishes 
as z ^ 0. 

In Figure 4 we plot Tc{z) and Tp{z) at fixed time delay, il- 
lustrating the two dynamical regimes that we observe. In the 
limit that Cc « Cp, the 2TM predicts that Tc{z) = Tp{z) 
(blue line) throughout the film. The temperature of both fluids 
varies on a length scale determined by the film thickness. In 
the opposite limit, Cc » Cp, which we believe is applicable 
to SRO, Tc{z) and Tp{z) are equal and nearly constant except 
near the interface, where Tc remains constant while Tp ap- 
proaches T. In both regimes, ATi decays exponentially with 
time, although the scaling of r with L is different in two cases, 
as we explain more fully below. 

To understand these simulations, we note that the character- 
istic time for the two-fluids to reach quasiequilibrium is given 
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FIG. 4: Temperature profiles for confined fluid and phonons for (a) 
Cc ^ Cp and < Kp- (b) Cc ^ Cp and Kc > t^P- The tempera- 
ture profile near the surface {z — 0) is dominated by the phonons in 
(a), and by confined fluid in (b). The shapes of these profiles remain 
constant throughout the decay. In both cases, Hip in the substrate is 
greater than that in the film. 



by Teq = Ceff/G, where Cg// = CcCp/C. On this time- 
scale, thermal energy stored in the phonons propagates a dis- 
tance Lp = ^{CclCfl^. The length scale J = (A^p/G)^/^ is 
the disequilibrium length identified previously |[lllll2|]. When 
Cc « Cp, we have Lp « S, and the temperatures of the 
two-fluids are locked together throughout the film. In the op- 
posite regime, where most of the energy is stored in the con- 
fined DOF, Lp S. 

While our simulations indicate that the decay of ATi with 
time is exponential in both regimes, we find different scal- 
ing with L; r (X when the phonons dominate the specific 
heat, whereas r ex L when the magnons dominate. The con- 
trasting spatial profile of Tp is the key to the different scaling. 
When Cc « Cp, Tp{z), r ^ L'^/Dp, where D is the ther- 
mal diffusivity. In the opposite limit the phonon heat current 
J is approximately tZpTp{0)/S (where hZp is the phonon ther- 
mal conductivity) within a length S from the interface and is 
essentially zero elsewhere. The thermal time constant is gov- 
erned by the relation r = Q/Q, where Q is thermal energy 
per unit area. Substituting Q = J and Q = CcLTp{0) we 
obtain. 
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(3) 



Eq. 4 is the key result that emerges from the 2TM model 
in the Cc » Cp regime. While the linear dependence of r 
on L indicates the existence of TBR, Eq. 4 is derived from a 
model that assumes perfect phonon transmission at the inter- 
face. The effective TBR seen here is the same as predicted to 
arise from disequilibrium between free and confined DOF's, 
for steady-state heat flow liL J2]. It is somewhat surprising 
to see the same effect in time-dependent heat flow, where it 
has been argued that two-fluid disequilibrium does not influ- 
ence the rate of cooling when t > Teq. However, our analysis 



shows that the presence of the interface leads to persistence of 
disequilibrium on the much longer time scale, Teq{L/6). 

Next, we consider the dependence of r on L and T in 
the light of the discussion above. To compare the mea- 
sured t{L) with Eq. 4, we assume that the effective phonon 
mean-free path, lp is given approximately by the relation 
= l~Q + where Ipo is the mean-free path in the 
absence of interfaces. Substituting the effective thermal con- 
ductivity K,plp/lpo into Eq. 4 yields r{L) = a{L -\- Ipo) 
for L >> Ipo and r(L) = a{lpoLy^'^ for L << Ipo, 
where a = {Cc/GtipY^'^ . The predicted linear to square 
root crossover in the dependence of r on decreasing L is fully 
consistent with the experimental data presented in Fig. 3b. 
We identify the crossover value of L of approximately 15 nm 
with the phonon mean-free path at 5 K. To make further quan- 
titative comparison with experiment, we note that from the 
measured value a=-2.1 ps/nm and literature value C = 10^ 
W/m-K, we obtain = 4 x 10^ W/m^-K. To obtain G 

we need K,p, which is not easy to extract from experiments, 
as they measure the total hz. Nevertheless, we can base an 
estimate on the kinetic theory relation, hZp = CpVslp\ substi- 



tuting /p = 15 nm. 



10^ m/s, and C„ = 24 W/m^-K 



(based on 6d =370 K Q) yields K.p ^1.2x10"^ W/m-K. 
This estimate, in turn, yields G 1.6 x 10^^ W/m^-K, which 
is consistent with equilibration rates seen in other magnetic 
materials I9i[il[ii]. 

Finally, we return to our original motivation for consider- 
ing disequilibrium effects, the weak temperature dependence 
of r for thick films. Above we have shown that in the thick 
film limit r/L ^ {Cc/Cy\Gvslp)-^^'^ . According to this 
formula, there is an accidental cancelation of the strongly T- 
dependent terms; the T^/^ dependence of the magnon spe- 
cific heat cancels the square root of the dependence of the 
phonon specific heat. The residual weak T dependence of 
t/L (see Fig. 3c) is associated with the term {Gvslp)~^^'^ . 
For thin films, the increase in lpo{T) with decreasing T is cut- 
off at L and below this T, r/L becomes simply proportional 
to G~^l^ . Again referring to Fig. 3c, and focusing on the 5 
nm sample, we associate the change in slope at 50K with the 
crossover between /^o > ^ and /^o < ^ regimes. 

To summarize, we have used TRMOKE to time-resolve the 
local T in thin films of the ferromagnet SRO. The scaling of 
Ti with laser intensity and substrate T shows that the specific 
heat of SRO is dominated by the magnon contribution in a 
broad T range. The flow of heat from film to substrate ap- 
pears to be limited by an effective thermal boundary resistance 
arising from disequilibrium between confined and propagat- 
ing degrees of freedom. 
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